<?php
/**
 * 国际坐标和国内坐标转换处理
 * WGS84坐标系，即国际坐标系由美国定义（我国的北斗也是采用此坐标协议），谷歌地图采用
 * GCJ02坐标系，即火星坐标系，由我国国测局制定的标准，是对WGS84的加密，高德地图、腾讯搜搜地图、阿里云地图、灵图51ditu地图、谷歌中国地图
 * BG-09坐标系。即百度坐标系，GCJ02坐标系经加密后的坐标系
 */

namespace service;

class Gps
{
    const x_PI = 52.35987755982988;         // π * 3000.0 / 180.0  这里π = 3.14159265358979324
    const PI   = 3.1415926535897932384626;  // π
    const a    = 6378245.0;                 // WGS长轴半径
    const ee   = 0.00669342162296594323;    // WGS偏心率的平方

    /**
     * WGS84转GCj02(北斗转高德)
     * @param int|string $lng   经度
     * @param int|string $lat   纬度
     * @return array
     */
    static public function beiDouToAMap($lng = 0, $lat = 0)
    {
        if (self::outOfChina($lng, $lat)) {
            return array($lng, $lat);
        } else {
            $dlat      = self::transFormLat($lng - 105.0, $lat - 35.0);
            $dlng      = self::transFormLng($lng - 105.0, $lat - 35.0);
            $radlat    = $lat / 180.0 * self::PI;
            $magic     = sin($radlat);
            $magic     = 1 - self::ee * $magic * $magic;
            $sqrtmagic = sqrt($magic);
            $dlat      = ($dlat * 180.0) / ((self::a * (1 - self::ee)) / ($magic * $sqrtmagic) * self::PI);
            $dlng      = ($dlng * 180.0) / (self::a / $sqrtmagic * cos($radlat) * self::PI);
            $mglat     = $lat + $dlat;
            $mglng     = $lng + $dlng;

            return [$mglng, $mglat, 'lng' => $mglng, 'lat' => $mglat];
        }
    }

    /**
     * GCJ02 转换为 WGS84 (高德转北斗)
     * @param int|string $lng   经度
     * @param int|string $lat   纬度
     * @return array
     */
    static public function aMapToBeiDou($lng = 0, $lat = 0)
    {
        if (self::outOfChina($lng, $lat)) {
            return array($lng, $lat);
        } else {
            $dlat      = self::transFormLat($lng - 105.0, $lat - 35.0);
            $dlng      = self::transFormLng($lng - 105.0, $lat - 35.0);
            $radlat    = $lat / 180.0 * self::PI;
            $magic     = sin($radlat);
            $magic     = 1 - self::ee * $magic * $magic;
            $sqrtmagic = sqrt($magic);
            $dlat      = ($dlat * 180.0) / ((self::a * (1 - self::ee)) / ($magic * $sqrtmagic) * self::PI);
            $dlng      = ($dlng * 180.0) / (self::a / $sqrtmagic * cos($radlat) * self::PI);
            $mglat     = $lat + $dlat;
            $mglng     = $lng + $dlng;

            return array($lng * 2 - $mglng, $lat * 2 - $mglat);
        }
    }

    /**
     * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换 百度转高德
     * 即 百度 转 谷歌、高德
     * @param string $bd_lon    百度经度
     * @param string $bd_lat    百度纬度
     * @return array
     */
    static public function baiDuToAMap($bd_lon = '', $bd_lat = '')
    {
        $x      = $bd_lon - 0.0065;
        $y      = $bd_lat - 0.006;
        $z      = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * self::x_PI);
        $theta  = atan2($y, $x) - 0.000003 * cos($x * self::x_PI);
        $gg_lng = $z * cos($theta);
        $gg_lat = $z * sin($theta);

        return array($gg_lng, $gg_lat);
    }

    /**
     * GCJ-02 转换为 BD-09  （火星坐标系转百度 即谷歌、高德转百度）
     * @param string $lng   经度
     * @param string $lat   纬度
     * @return array
     */
    static public function aMapToBaiDu($lng = '', $lat = '')
    {
        $z      = sqrt($lng * $lng + $lat * $lat) + 0.00002 * sin($lat * self::x_PI);
        $theta  = atan2($lat, $lng) + 0.000003 * cos($lng * self::x_PI);
        $bd_lng = $z * cos($theta) + 0.0065;
        $bd_lat = $z * sin($theta) + 0.006;

        return array($bd_lng, $bd_lat);
    }

    static private function transFormLat($lng = '', $lat = '')
    {
        $ret  = -100.0 + 2.0 * $lng + 3.0 * $lat + 0.2 * $lat * $lat + 0.1 * $lng * $lat + 0.2 * sqrt(abs($lng));
        $ret += (20.0 * sin(6.0 * $lng * self::PI) + 20.0 * sin(2.0 * $lng * self::PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($lat * self::PI) + 40.0 * sin($lat / 3.0 * self::PI)) * 2.0 / 3.0;
        $ret += (160.0 * sin($lat / 12.0 * self::PI) + 320 * sin($lat * self::PI / 30.0)) * 2.0 / 3.0;

        return $ret;
    }

    static private function transFormLng($lng = 0, $lat = 0)
    {
        $ret  = 300.0 + $lng + 2.0 * $lat + 0.1 * $lng * $lng + 0.1 * $lng * $lat + 0.1 * sqrt(abs($lng));
        $ret += (20.0 * sin(6.0 * $lng * self::PI) + 20.0 * sin(2.0 * $lng * self::PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($lng * self::PI) + 40.0 * sin($lng / 3.0 * self::PI)) * 2.0 / 3.0;
        $ret += (150.0 * sin($lng / 12.0 * self::PI) + 300.0 * sin($lng / 30.0 * self::PI)) * 2.0 / 3.0;

        return $ret;
    }

    /**
     * 判断是否在国内，不在国内不做偏移
     * @param int $lng  经度
     * @param int $lat  纬度
     * @return bool
     */
    static private function outOfChina($lng = 0, $lat = 0)
    {
        return ($lng < 72.004 || $lng > 137.8347) || (($lat < 0.8293 || $lat > 55.8271) || false);
    }

}